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Using a bit-string model of evolution, we find a successful route to diploidy and sex in simple 
organisms. Allowing the sexually reproducing diploid individuals to also perform mitosis, as they 
do in a haploid-diploid cycle, leads to the complete takeover of the population by sexual diploids. 
This mechanism is so robust, that even the accidental conversion and pairing of only two diploids 
give rise to a sexual population. 
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I. INTRODUCTION 

The evolution of sex through Darwinian selection in 
spite of the seeming odds with respect to asexual re- 
production through cloning, is a still unresolved prob- 
lem which has attracted a great deal of recent 
activity [ . We would like to present here results from 
a very simple model based on unicellular organisms where 
sex seems to have first emerged, without any complicat- 
ing factors such as the cost of male-female differentiation, 
parental care, mating preferences, maternity periods and 
so on, which we assume only to have arisen in higher 
organisms. These factors may be responsible for the dif- 
ferentiation, modification or even at times reversal (as 
in meiotic parthenogenesis |^,p.234) of the patterns that 
may have become accidentally established in the remote 
evolutionary past. Unicellular organisms are not only 
simpler to model, but this choice also imposes more strin- 
gent conditions on the establishment of a sexual popula- 
tion, since it severly limits the number of viable offspring 
that a pair of parent cells can have. By confining our- 
selves to "worst case scenarios" we hope to be able to 
obtain plausible lower bounds to the feasibility of sex as 
a reproductive mechanism. 

Sexual reproduction as we know it today in higher 
plants and animals can be regarded as a haploid-diploid 
cycle with a highly abbreviated haploid phase, where the 
(haploid) gametes typically do not perform mitosis. ||] 
These haploids are not viable unless two of them fuse and 
once more make a diploid cell, which then multiplies by 
mitosis and eventually forms the mature individual. 

The haploid diploid cycle (HDC) where, under un- 
favorable conditions the population becomes diploid, is 
found in many species ^, and was the motivation 
behind the hypothesis of Jan, Stauffer and Moseley ||] 
(JSM) who proposed that diploidy and sex may have 
first arisen as a way to escape death, when a simple, uni- 
cellular individual is threatened by too many deleterious 
mutations. The fitness is taken to be a step function 
of the number of mutations, also setting the threshold 
for conversion to diploidy and sex. Here, the number of 
deleterious mutations can be read as the distance from 



the ideal case, called the wildtype, which can be altered 
as a result of a change in the environment, so that the 
individual is not as well adapted as before the change. 

In previous work ||, |l^ we showed that the JSM hy- 
pothesis indeed leads to a steady source of diploid sexu- 
als, given a population of haploids which multiply by mi- 
tosis. This leads to a steady state distribution of coexist- 
ing haploids and diploids in a constant population. The 
premises adopted in [0 were very restrictive. Two sex- 
ual "parent" cells were allowed to have only one offspring, 
after which they died. In case one considers greater num- 
ber of offspring, by allowing a greater number of viable 
gametes to be formed, as is indeed possible in many uni- 
cellular organisms one finds that the diploid, sexual 
population completely takes over. 

In the same paper we have also considered the sit- 
uation where conversion to sex occurred with a constant 
probability a (which we varied between and 1), inde- 
pendently of the fitness of the individual. This strategy 
also gave rise to a sexual population making up a small 
but macroscopic fraction of the total, so that the JSM 
threshold mechanism for conversion to sex did not prove 
to be necessary, although it was more successful. On the 
other hand the JSM hypothesis has the additional at- 
traction of providing a mechanism which could trigger 
the fusion of two haploid cells to form a diploid: it is 
known 0|,p.l49, that extensive genetic damage can lead 
to gene repair via genetic transmission between two hap- 
loid cells; the fusion of two haploids could be seen as an 
extreme form of such behaviour. An alternative means 
of forming a diploid from a haploid cell is via endomi- 
tosis ||],p.230, as in the first step in meiotic partheno- 
genesis. Endomitosis is the process whereby the genetic 
material in a cell is duplicated without subsequent cell di- 
vision as in normal mitosis. Again, it is not implausible 
that this process originated as a result of grave genetic 
damage which precluded the successful completion of mi- 
tosis. 0, p. 149. 

In the present paper we show that, if the diploid cells, 
once formed, are also allowed to multiply by mitosis, as 
indeed they do in a haploid-diploid cycle, the whole pop- 
ulation is taken over by diploid cells, who perform facul- 
tative sex if they are once more threatened by extinction 
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due to too many deleterious mutations. Moreover we 
show that even an episodic conversion to sex, involving 
as few as only two individuals who survive to mate, leads 
to a steady state made up solely of sexual types. 

In the next section we briefly describe our algorithms 
and report the simulation results. In section 3, we pro- 
vide a discussion and some pointers for future research. 



II. BIT-STRING MODEL FOR THE 
CONVERSION TO SEX - ALGORITHMS AND 
SIMULATION RESULTS 

The type of model we consider here is the same as in 
Ref. ||l^. Below and in Section A, we recall the defini- 
tions. In Section B and C wc introduce new rules to test 
the autonomous viability of sexual populations. 

Each haploid one-celled organism consists, for our pur- 
poses, of a 15-bit string of "0"s and "l"s, representing the 
genetic code in a 16-bit computational word. |l^, |l2|, |l^ 
We use the bit defining the "sign", to specify whether 
the individual is asexual or sexual (-). A muta- 
tion consists of flipping a randomly chosen bit except 
the sign bit, with a constant probability, F, for each 
individual per generation. Since the genetic difference 
between individuals of the same species is typically less 
than 10% [|l|,p.52, this rather short string for the genetic 
code may be considered as a coarse grained model for the 
complete genome of the individual, which we divide up 
into different zones, retaining a "0" where there are no 
mutations, and flipping the bit to "1" if there are one or 
more mutations in this zone. The wild type is a string of 
all "0"s. Therefore, at each locus, a "1" corresponds to a 
deleterious mutation (which we will call "mutation," for 
short, where this is not liable to lead to any confusion.) 
Diploid organisms have two bit strings, which are allowed 
to be different. The number of deleterious mutations m is 
simply the number of "l"s for a haploid individual. For a 
diploid, the number of "expressed" deleterious mutations 
is the number of loci at which both homologous alleles 
are set to "1," i.e., we assume that deleterious mutations 
are recessive. 

The total population is fixed at iV = 10'^, and we have 
chosen F = 1/A^. This corresponds to a mutation rate 
per allele per generation of 6 x 10~^, which is com- 
parable to the typical mutation rates encountered in eu- 
karyotes |H,p.64. 

The probability of survival (or fitness) as a function of 
m is given by a Fermi-like distribution |jlj|, P{m), 



P{m) 



exp[/3(m - /x)] + 1 



(1) 



For large (3 (or "low temperatures," in the language of 
statistical mechanics), P(m) behaves like a step function. 
Individuals with m > fi die, those with m < ^ survive, 
and those with m ~ /i survive with a probability of 1/2. 
In the simulations we confined ourselves to low temper- 
atures (/3 = 10). We chose fj, = 4 which just allows us 



enough variability without leading to totally unrealistic 
mutational loads. 

We start with a set of N asexual (haploid) wildtypes. 
In each generation about FA^ individuals suffer muta- 
tions; they are killed off or retained according to the fit- 
ness function (|^), and the population restored to N by 
once cloning as many randomly chosen survivors as nec- 
essary. The population of haploid asexuals settles down 
to a minimally stable steady state distribution (see 
Table I) as shown in Fig. 1, independently of the value 
of F, for F > 1/N. g ^ 

For comparison, we have also performed simulations on 
a diploid population reproducing asexually, according to 
the same rules as stated in the previous paragraph. We 
have found that the diploid population reaches a steady 
state with an m distribution peaked at m — 2, rather 
than 771 = 3 as found for haploids. (See Fig. 1 and Table 
I)- 



Conversion to sex 



Let us briefly summarise our algorithm for conversion 
to sex, which we also used in our previous study (see 
Model A, Ref. |o|). 

We choose the steady state of the haploid population 
as our initial state from which to start the conversion 
to sex, as that would be a most likely "natural state" 
encountered at this stage. Once the haploid population 
reaches a steady state, we allow those individuals that 
are threatened by extinction to convert to diploidy and 
sex, by implementing the following rule at each genera- 
tion: If an asexual individual with m = 4 has survived, 
it is converted to an active sexual, by deterministically 
and irreversibly switching its sign bit to (-). It performs 
endomitosis and becomes diploid. If there are already 
sexual, diploid organisms in the population, they will 
also be made active if to, the number of their expressed 
mutations, is > 4 (otherwise they do not participate in 
the reproduction cycle; hence the conversion to sex may 
be termed non- hereditary) . Finally, all the active sex- 
ual organisms pair randomly and engage in sexual repro- 
duction, where they each contribute one gamete (formed 
via one step meiosis) towards a single diploid sexual off- 
spring. If we denote the genotypes of parents as {Aa} 
and {Bb} respectively, then the genotype of the offspring 
is cither {AB}, {Ab}, {aB} or {ab}. No crossover occurs 
during this one step meiosis. If there is only one active 
sexual at a certain time step then it must wait subsequent 
generations until it either finds a partner or it dies. We 
keep the population constant [p^ by cloning randomly 
selected haploid individuals to make up the deficit. 

(The consequences of the conversion to hereditary and 
obligatory sex can be found in Ref. pC|].) 



3 



B. Mitotic Diploid Sexuals Win over Haploids 



Episodic conversion to diploidy and sex leads to 
sexual steady state 



In the previous subsection we described a scenario 
where too big a mutation load meant conversion to 
diploidy and sex, if the individual survived. It should 
be realised that in this scheme the haploids, who mul- 
tiply by cloning, provide a steady source for the diploid 
sexuals, whose numbers are halved every time they mate. 
Now we want to test whether the diploid, sexual popu- 
lation can survive autonomously, if the sexuals are also 
allowed to perform mitosis themselves. 

Once a steady state with coexisting haploid and diploid 
populations is reached via the algorithm described above, 
we switch off the conversion of haploid organisms to 
diploidy and sex. The diploid individuals mate when 
they face extinction due to too many mutations. We 
now keep the total population constant by making up 
for the deficit in the population at each step by cloning 
randomly selected individuals, regardless of whether they 
are sexual or asexual. In colonies undergoing a haploid- 
diploid cycle, it is quite frequently the case that the 
diploid phase of the cycle also involves multiplication by 
mitosis. 

The result is that the diploid, sexual individuals com- 
pletely win over the population. The haploids which now 
can not compete with the sexuals, become extinct. Here 
we find that the haploid phase of the HDC becomes ab- 
breviated to the point where haploids appear only as ga- 
metes which do not perform mitosis. This is exactly the 
situation in highly evolved sexual organisms. 
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FIG. 1: The steady state distribution of pure haploids with 
respect to the number deleterious mutations m, compared 
with the distribution of asexual diploids over the number of 
expressed mutations m. See text. The histograms are nor- 
malized to unity. 



One can also think of a scenario in which conversion to 
sex takes place accidentally over a short period of time af- 
ter the asexual population settles down to a steady state. 
In fact this is probably the most realistic situation, given 
the random nature of the mutations. 

The way we actually implemented this in the computer 
code was by deterministically switching the sign bits of 
the first two asexual individuals to survive with m = 4, 
and then turning off the possibility of further conversion. 
These then form two sexual individuals by endomitosis, 
and if the first survives long enough so that it can mate 
with the second, will give rise to one sexual offspring. 

The rest of the rules are as explained before; in each 
generation we clone randomly chosen individuals to make 
up the deficit in the population, regardless of whether 
they are asexuals or sexuals. We allow the diploid indi- 
viduals to mate when they face extinction in the course 
of their lives. 

Suprisingly, the sexuals capture the population in 95 
percent of the performed runs. (In the rest, the single 
diploid, which is still at the threshold of extinction with 
an m value of 4, may not survive until a partner arrives.) 



D. Relative Fitness 

The steady state m distributions for the sexual pop- 
ulations treated in Section II. B and C are much better 
than the haploid distribution. However, they are iden- 
tical with the distribution for asexual diploids. (see Ta- 
ble I and Fig. 1) This means that as soon as diploidy 
is achieved, it is so successful in screening the effects 
of deleterious mutations that once diploid, the organ- 
ism practices sex very infrequently, making these three 
distributions identical. 

The diploid m-distribution in Fig. 1 should be con- 
trasted with the results of Ref. [Q (Model A, non- 
hereditary and hereditary; also see Subsection A, above) 
where cloning of the sexuals is not allowed, but there is a 
steady influx of new sexuals from the haploids (see Table 
I). In the non-hereditary model, the m distribution has 
shifted nearer the wildtype than in the purely diploid case 
(although the peak is again at m = 2); in the hereditary 
model, the peak has shifted to to — 1. Thus, the fre- 
quency with which sex is practiced has a salutary effect 
on the distribution of expressed deleterious mutations in 
the population. 

On the other hand it should be noted that, with the 
Fermi-like fitness function (Eq.(l)) for large f3, to which 
we have confined ourselves, the survival probability is ~ 
constant for m < /i, and therefore does not discriminate 
between steady state populations that differ solely in the 
shape of their m distributions for m < ii. A better mea- 
sure might be the deathrate, namely, the average number 
of individuals which are eliminated at each step, under 
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TABLE I: The distribution of expressed deleterious mutations and death rates for the models discussed in Section II. The 
m-distributions are percentages (±1%) computed over each type if coexisting; the death rates are per capita per generation. 
(* see Ref. @) 
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the criterion given in Eq.(l). 

Our findings, which are rather revealing, are summa- 
rized in Table I. Since the fluctuations were very large, 
we have only reported order of magnitude results for the 
death rates. Pure diploids have the lowest death rate, an 
order of magnitude smaller than pure haploids. Mitotic 
(including the "episodic" case) sexuals are very close to 
diploids. For non-hereditary sexuals the rates are com- 
parable to asexuals, but hereditary sexuals have a death 
rate that is two orders of magnitude larger than asexuals. 
We infer from this that the success of hereditary sexuals 
in lowering their mutation load is due to the much greater 
rate at which they can select-out highly mutated genes 
through death. 

Thus it is seen that the clear advantage of sex over asex 
(haploid or diploid) in our model can only become mani- 
fest with a sufficiently high frequency of sexual reproduc- 
tion, which is driven by the mutation rate F, and more- 
over in situations which can discriminate between the 
different adaptabilities (smaller typical m) of the differ- 
ent modes of reproduction, such as time varying environ- 
ments. jl7| This is similar to the finding of Pekalski [^ , 
who has considered environments (wildtypes) that vary 
over time, and found that the benefits of sexual repro- 
duction are enhanced by higher birthrates. For lower 
birthrates, meiotic parthenogenesis, which can be com- 
pared to diploid asex, has a slight advantage over sex, a 
difference from our results which can be ascribed to the 
exponential fitness function he uses. 



III. DISCUSSION 

This paper is a culmination of a series of studies where 
we have considered very stringent rules for the survival 
and propagation of diploid, sexual individuals in com- 
petition with haploids. Here we have finally been able 
to show that a pair of simple, unicellular organisms who 
have accidentally converted to diploidy, and which sub- 
sequently engage in sexual reproduction, begetting one 
sexual offspring, can give rise to a population of sexual 
types which totally take over a finite population, pro- 
vided they are also allowed to multiply by mitosis, on an 



equal footing with the haploids in the population. They 
engage in sex when the going gets tough, that is, when the 
number of their expressed deleterious mutations exceeds 
a certain number. This success seems to vindicate the 
hypothesis of JSM that sex could have been a mech- 
anism of last resort when simple organisms were faced 
with extinction. 

It is interesting to pose the question of what happens 
if we do not allow the diploids to clone themselves, but, 
on the other hand, allow the gametes, under special con- 
ditions, to enter a haploid phase where they multiply by 
mitosis. Eventually these haploids will be allowed once 
more to fuse and give rise to diploids. Within the present 
scheme, since the haploids convert to diploidy only at the 
threshold, m = 4, the gametes of the diploid individuals 
have TO > 4, i.e., they are not viable. Therefore the 
extension of the haploid phase calls for a modification of 
the rules in such a way as to allow the gametes to survive 
nevertheless, e.g., by forming "spores." 

One way the parameters of the present model could be 
modified is to take the threshold for survival to be differ- 
ent (greater than) the threshold for the conversion to (or 
practice of) sex. Then, if a haploid gamete were to be 
considered as a member of the haploid population, and 
subjected to the same rules for selection and reproduc- 
tion, it would seem as if one would get haploid phases of 
arbitrary duration. However, upon closer inspection we 
see that this is not true: Since the gamete comes from a 
diploid organism which engages in sex due to an excess 
mutational load, it already has deleterious mutations in 
excess of (or at least equal to) the sex threshold. (For the 
diploid organism to have m expressed mutations, both of 
its genetic strings must have at least that many). So, ac- 
cording to the rules which we have adopted, the abbrevi- 
ated haploid cycle is automatically selected, since within 
this scheme the gamete would immediately seek a mate. 
The lengthening of the haplophase would call for further 
evolution, e.g., such mutations as would lead to the en- 
hancement of the fitness of the haploid gametes with m 
in excess of the sexual threshold (say by the formation of 
spores, to survive in an environment which has become 
disfavorable [3,4]). Since this will provide a mechanism 
for the diploids once more to make up for the reduction 
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in their numbers due to the single-offspring mating rule, 
we expect this HDC to win out over the haploid, asexual 
population. 

To understand the seemingly contradictory of results 
from similar models investigating the phenomenon of the 
selection of sex as the dominant mode of reproduction, 
it is of interest to compare our model with that of Red- 
field ||l6| and Cui et al. , . The former employs a more 
complex system in which there is sexual differentiation 
between males and females, with the additional feature 
of a much greater mutation rate for males than for fe- 
males, eventually making meiotic parthenogenesis safer 
than sex. Redfield has used three different "selection 
functions" (survival probabilities), falling off exponen- 
tially (no "epistasis," or correlation between the effects of 
successive mutations) quadratically (some positive epis- 
tasis) or as a step function (extreme positive epistasis) 
with r7i[^. The fact that the step-function ("trunca- 
tion") survival probability alone leads to a much better 
fitness resulting from sex, tells us that the form of the 
survival probability which is adopted determines the out- 
come very strongly. The step function favors sex because 
it does not punish small deviations from the wildtype, 
thereby allowing genetic diversity, while it severely pe- 
nalizes large deviations, helping eliminate highly mutated 
individuals. 

This insight enables us to understand the results of 
Cui et al. who have considered a diploid population re- 
producing either asexually or sexually, depending upon 
a fixed probability. There is a constant mutation rate. 
The survival probability is chosen to be of the "indepen- 
dent mutation," i.e., exponential type. These authors 
find that the diploid population accumulates such a large 
number of deleterious mutations (which however do not 
get expressed), that sexual reproduction, i.e., the ran- 
dom pairing of gametes from different parents, results 
in a disasterous reduction in the fitness of the offspring. 
Moreover, sexually reproducing diploid populations are 
susceptible to invasion by asexually reproducing ones, 
while the converse is not true, i.e., asexually reproducing 
diploid populations are not invaded by sexually repro- 
ducing ones. These results change drastically in favor of 
sex when they introduce cell senescence, which means a 
cell can clone itself only a certain number of times be- 
fore it stops dividing and dies, unless it engages in sexual 
reproduction, which resets the senescence clock. 

We would like to argue that in the model of Cui et 
al., introducing cell senescence is very much like turning 
on a "truncation like" survival function, which, as in our 
model, also controls the switch to sexual reproduction. 
In the presence of a constant mutation rate, the num- 
ber of times a cell has cloned is another way to measure 
its (mean) mutation load, which, in an asexual haploid 



population of course gets directly expressed, whereas for 
diploid individuals hides behind the dominant unmutated 
alleles. Introducing a cell senescence threshold, beyond 
which the individual is either killed off or has to engage in 
sexual reproduction, is therefore equivalent, in an average 
way, to the JSM criterion, applied to a diploid popula- 
tion. Seen in this way, it is very gratifying that the Cui 
et al. model in fact corraborates our findings regarding 
the greater fitness of diploid populations that frequently 
engage in sex, given a step-function like survival proba- 
bility. 

Very recently, two variants of the Redfield model, 
which assign "harsher selection" to males (rather than 
relatively higher mutation rates) have been consid- 
ered ^ . These studies both find that the mutational 
load on the whole population decreases as a result, and 
that the relative fitness of sexual females is increased by 
more than enough to compensate for the twofold cost of 
sex in anisogamous populations, at least for extremely 
high values of the average mutation rate. Agarwal finds 
that this can happen even in the absence of synergystic 
epistasis, in his model. Wilke et. al j23j have meanwhile 
found that high mutation rates select for low replication 
rates and flat regions of the fitness surface. 

It should be noted that the Fermi-type fitness func- 
tion (1) adopted here|^, ^, extrapolates, for finite 
temperatures (smaller values of /3) between highly syner- 
gistic (step function) and independent (exponential) sur- 
vival probabilities 124] as a function of m. A recent study 
by Peck and Waxman indicates that competition for 
limited resources can lead to synergy between successive 
mutations, leading to truncation, or step function-like 
survival probabilities, which they also find favors sexual 
populations. 

In this paper we have presented numerical results for 
a model involving unicellular organisms, which may shed 
light on how sex and diploidy emerged and established a 
foothold in the protozoon world. Many different mecha- 
nisms have so far been proposed via which sex may prove 
advantageous or otherwise, in more highly evolved organ- 
isms. We believe that care should be taken while propos- 
ing any single mechanism, such as an adaptation to resist 
infestation by parasites (26[ , for the preferance of sex 
by organisms that range from the unicellular eukaryotes 
to trees or human beings psf . More complex organisms 
may have elaborated much more complex survival mech- 
anisms and behaviorial patterns, which stabilize or desta- 
bilize already evolved traits. 
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